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ABSTRACT 

We present a census of molecular outflows across four active regions of star formation in 
the Perseus molecular cloud (NGC 1333, IC348/HH211, L1448 and L1455), totalling an area 
of over 1000 arcmin 2 . This is one of the largest surveys of outflow evolution in a single molec- 
ular cloud published to date. We analyse large-scale, sensitive CO / = 3 — > 2 datasets from 
the James Clerk Maxwell Telescope, including new data towards NGC 1333. Where possible 
we make use of our complementary 13 CO and C 18 data to correct for the 12 CO optical depth 
and measure ambient cloud properties. Of the 65 submillimetre cores in our fields, we detect 
outflows towards 45. 24 of these are marginal detections where the outflow's shape is unclear 
or could be confused with the other outflows. We compare various parameters between the 
outflows from Class and I protostars, including their mass, momentum, energy and momen- 
tum flux. Class outflows are longer, faster, more massive and have more energy than Class 
I outflows. The dynamical time-scales we derive from these outflows are uncorrected to the 
age of the outflow drivi ng source, com puted from the protostar's bolometric temperature. We 
confirm the results of Bontemps et al.| that outflows decrease in force as they age. There is a 
decrease in momentum flux from the Class to I stage: (Fqo) = (8 ± 3) x 10 -5 compared to 
(1.1 ±0.3) x 10 -5 M kms -1 yr _1 , suggesting a decline in the mass accretion rate assuming 
the same entrainment fraction for both classes of outflow. If F ra( j = L^x/c is the flux expected 
in radiation from the central source, then fco(ClassI) ~ lOO/^ad and Fqo (Class 0) ~ 1000F ra d. 
Furthermore, we confirm there are additional sources of mass loss from protostars. If a core's 
mass is only lost from outflows at the current rate, cores would endure a few million years, 
much longer than current estimates for the duration of the protostellar stage. Finally, we note 
that the total energy contained in outflows in NGC 1333, L1448 and L1455 is greater than 
the estimated turbulent energy in the respective regions, which may have implications for the 
regions' evolution. 

Key words: submillimetre - stars: formation - stars: mass-loss - ISM: jets and outflows - 
ISM: individual: Perseus. 



1 INTRODUCTION 

Molecular outflows are ubiquitous in star-forming regions, with ev- 
ery protostar thought to undergo a period of outflow activity (e.g. 
|Shu, Adams & L izano 1987 ). Outflow observations provide unique 
information about the mass-loss history of the parent protostar and 
thereby the accretion process which determines its evolution. Re- 
cently there has been renewed interest in the suggestion that out- 
flows may drive turbulence in molecular clouds (Norman & Silk| 
1980 ), given that the total energy they contain may be a significant 
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portion of a molecular cloud's turbulent energy and they can cause 



much disruption of their natal environment (e.g. |Knee & Sandell 
|2000[ | Wolf-Chase, Barsony & O , Linger|2000||Stanke& Williams 



2007 ). Theoreticians are divided about whether an efficient mech- 
anism exists to feed energy from a small quantity of fast-moving 
material in an outflowing jet into random isotropic motions on 



molecular cloud scales (e.g. Cunningham et al. 2006 2009 versus 
Banerjee, Klessen & Fendt 2007}. However, analyses of data from 
regions of molecular clouds with high densities of outflows sug- 
gest that outflows are not the dominant driving mechanism ( [Brunt, | 
|Heyer & Mac Low|[2009l |Padoan et~aT]|2009] l. If they were, then 
cloud lifetimes should be short (consistent with the rapid forma- 
tion scenario proposed by Hartmann, Ballesteros-Paredes & Bergin 
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200 1|> or the s tar formation rate very high (see Williams, Plam-| 



beck & He yer 2003). In addition, the detection of outflows from 



candidate star-forming cores, located through say (sub)millimetre 
continuum surveys (e.g. Ward-Thomps on et al.| |2007), offers a 
way to classify such sources as either starless or protostellar (e.g. 
|Hatchell, Fuller & Richer|2007b| hereafter [hFr| , which comple- 
ments methods based on their spectral energy distributions (SEDs, 
e.g. |Hatchell et al.|2007a||Enoch et al.|2009l >. 

Molecular outflows have diverse structures, reflecting both 
their driving sources and environments. Such diversity has lead 
to a wide range of plausible models (see |Arce et al.||2007| for a 
summary) with no single model accounting for all the features of 
outflows observed to date. Young, collimated outflows point to jet- 
driven models (e.g. |Raga & Cabrit|1993||Masson & Chernin|1993] >, 
yet older, wider-angle outflows evoke wide-angle wind models (e.g. 
|Shu et al.|1991| >. Some observational studies even suggest a two- 
component wind is necessary (e.g. Yu, Billawala & Bally 1999). 

In this paper we present a survey for molecular outflows to- 
wards the Perseus molecular clou d (hereafter simply Perseus) us 



Table 1. Details of regions included in the survey. The data towards IC348, 
L1448 and L1455 have been published previously ( Paper I). 



2 datasets |Curtis, Richer & Buckle|2010 



ing our 1 2 CO 7 = 3 
hereafter Paper I) and data presented here for the first time towards 
NGC 1333, the most active regions of star formation in Perseus. 
This outflow survey towards a single cloud is one of the largest to 
date, allowing us to look for statistical variations in outflow prop- 
erties. This work is the second in a series of papers surveying the 
kinematics of molecular gas in Perseus. Paper I presented the data 
acquisition and large-scale kinematic properties, whilst |Paper IIl| 
(Curtis & Richer 2010) expl ores the detailed gas motions inside 
star-forming cores. This study is arranged as follows: Section 2 
describes the observations we analyse, in particular the new data 
towards NGC 1333, and presents an overview of the survey before 
listing the naming conventions we apply to star-forming cores in 
our fields. In Section [3] we look for an outflow detection towards 
the position of every submillimetre core identified with SCUBA 
(Hatc hell et al.[ 2005 2007a) and analyse the corresponding detec- 
tion rates and core classifications. We set out the methods we have 
chosen to analyse the detected outflows in Section [4] (the outflows 
themselves are described in Appendi x[A| before exploring the evo- 
lution of these parameters in Section [5] Finally, we summarize our 
results in Section[6] 



2 OBSERVATIONS 

The 12 CO 7 = 3^2 datcQ we analyse are towards the four largest 
clusters of star-forming cores in Perseus (see Table[T]and Fig.[T]for 
an overview of the data). The datacubes of IC348/HH211, L1448 
and L1455 have been published previously ( Paper I) where we refer 
the reader for details of the reduction. Briefly, the data were taken 
as 'basket- weaved' scan maps using HARP (Heterodyne Array Re- 
ceiver Project; Buc kle et al.| [2009) on the James Clerk Maxwell 
Telescope (JCMT) over nine nights between December 2007 and 
January 2008. The data were taken in good to excellent observing 
conditions with a median system temperature of T sys = 325 K. The 
final data products are sampled on a 3 arcsec grid and are smoothed 
slightly spatially and spectrally to an equivalent full-width half 
maximum (FWHM) beam size of 16.8 arcsec (0.020 pc at 250 pc, 
our assumed distance to Perseus) and velocity channel width of 



Region 


Centre a 


Area 


G b 

u rms 






(h m s) 


(° ' ") 


(arcmin 2 ) 


(K) 


(kms" 1 ) 


NGC 1333 


03:28:56 


+31:17:30 


612 


0.07 


7.8 


IC348 


03:44:14 


+31:49:49 


199 


0.09 


8.9 


L1448 


03:25:30 


+30:43:45 


111 


0.04 


4.3 


L1455 


03:27:27 


+30:14:30 


127 


0.07 


4.5 



a Position of the map centre in J2000 coordinates. 

b Median rms noise across the map measured on spectra with lkms -1 
channel widths. 

c Characteristic velocity in the region, taken as the median value of the 
C ls O./ = 3— ^2 line centre velocity derived from Gaussian fits to every 
spectrum where the line is detected with a peak brightness temperature 
> 3^ (Paper I) . 



lkms -1 respectively. Smoothing to lkms -1 resolution provides 
additional sensitivity in the 12 CO line wings. 

Our survey area is more than doubled by including new 12 CO 
/ = 3 -+ 2 HARP data towards the young infrared (IR) cluster 
NGC 1333. In Fig. [2] we plot an integrated intensity image of 
the new data towards NGC 1333, which shows a wealth of out- 
flows criss-crossing the re gion. The largest previ ous CO / = 3 — > 2 



dataset was presented by 



Knee & Sandell 



(2000 



KSOOp 



hereafter 

and was only towards a small region (63 arcmin 2 ) in the cen- 
tre of our field. Our current dataset is nearly ten times larger in 
area (see Table [TJ. The data were taken over four nights between 
13th and 16th January 2007 and comprise around 27 hr of ob- 
serving time in total. The sky opacity recorded at 225 GHz var- 
ied from T225 = 0.025 to 0.125, with a median for all the scans of 
T 225 — 0.04, which contributed to the median system temperature 
of r sys = 210 K. The final mosaic is composed of multiple smaller 
tiles of various dimensions, which were chosen to cover not only 
the brightest 12 CO emission, but also all the 'green' nebulous fea- 
tures in the Spitzer IRAC map of the region. This explains why the 
map extends so far to the west. The idea was to make sure all of the 
outflow activity associated with NGC 1333 was covered in the sur- 
vey. Each tile is itself a basket- weaved scan map. The 4x4 detector 
array was inclined at 14deg to the scan direction as is customary 
to maximize the efficiency with HARP observations (Buckl e~et al.| 
2009). Unusually, the perpendicular spacing between adjacent scan 
rows was set to be 6 arcsec to yield Nyquist sampled observations 
(A/2D is 6.0 arcsec at 345 GHz) and ensure more of the 16 de- 
tectors cross the same spatial point in the sky, providing redun- 
dancy and smoothing-out variations in detector performance (i.e. 
noise). The backend correlator, ACSIS (Buc kle et al.| [2009), was 
configured to supply 1 GHz of bandwidth, split into 2048 spectral 
channels, each 488 kHz (0.42 km s -1 ) wide. The final NGC 1333 
datacube analysed is on a 6 arcsec spatial grid (distributed using a 
nearest-neighbour algorithm) and has been smoothed spectrally to 
a velocity resolution of 1 kms -1 as in the other regions. The re- 
sultant rms noise is listed in Table [T] All the standard observing 
procedures were followed at the telescope, in particular calibration 
spectra were taken frequently towards the CRL 618 standard to ver- 
ify its intensity was within the expected calibration tolerance. 

Any HARP data we present in figures is on the antenna tem- 



The transition frequency is 345.796 GHz. 



2 This was subsequently expanded by a number of surveys (Hatchell et al.| 
|2007b|[Hatchell & Dunham|2009) . 
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Figure 1. Overview of our outflow data across the Pe rseus molecular clou d. The central image of the entire cloud is C ls O / = 1 — >■ integrated intensity, 
scaled from to 4Kkms _1 from the FCRAO survey (Hatchell et al. 2005). The green contour encloses areas where the visual extinction is > 3 using data 
from the COMPLETE team (Ridge et al. 2006 ), whilst the orange contours signify the presence of high column densities of dust, marking where the SCUBA 
850 urn flux density is 200, 800 and 3200mJybeam- 1 jHatchell et al.|2005j . Our four survey regions (NGC 1333, IC348/HH211, L1448 and L1455) are 
marked with dark grey boxes. The enlargements of our four fields present 12 CO / = 3 — >■ 2 integrate d intensity, blu e- or red- shifted with respect to the ambient 
cloud velocity, contoured on SCUBA 850 \xm images, scaled from to 1 600 mJy beam -1 (Hatchell et al. 2005). The velocity limits of the integrations are 
the same as for Fig s . | A 1 1 |A7] (middle) , |A8 1 (top) and |A9| Except for IC348, the contours are at 5, 10, 20, 40, 80 and 160o\ where a is the rms noise on the 
integrated intensity image (computed per pixel). Towards IC348, where we have plotted lower-velocity outflows than the other regions (which trace the most 
spatial structure), the contours are at 3, 4, 5, 6, 7, 8, 9 and lOcr. 



perature scale (T£; Kutn er & Ulich|[T98T] ), which we convert to 
main beam brightness temperature, T m ^, where necessary in our 
calculations using T m ^ = T^/i]^. We use an efficiency, T7 m b = 
0.66, as measured during the commissioning of HARP. 



2.1 Survey overview and naming conventions 

Our 12 CO datacubes contain a large sample of outflows across 
Perseus, in total 23 Class and 14 Class I protostars are situated 
in our fields ( [Hatchell et aTl|2007a| ). The CO / = 3 2 line is 
arguably a better probe of outflows than lower / transitions (e.g. 
Takahashi et al. 2008) as it requires warmer temperatures for ex- 
citation (E v /k = 33.2 K above ground) that better match the gas 
temperatures observed in outflows (e.g. Hat chell, Fuller & Ladd| 
[T999l|Nisini et al.|2000) . Detailed images and a description of the 



outflows we find are presented in AppendixjA] In those descriptions 
and elsewhere we refer extensively to the submillimetre cores iden- 
tified by ( [Hatchell et al.|2005||2557a) from 850 \im SCUBA data. 
We label such cores with the prefix HFR, i.e. HFR1 is core 1 in 
|Hatchell et"aL] ( |2007a| ). These cores can be found in the SIMBAD 
database as [HRF2005]wz or [HRF2007]wi. A number of cores 
were originally identified in the Perseus Bolocam survey (Enoch 
|et al.|2006| and we refer to them using their original numbers with 
the prefix Bolo, e.g. BololOl. 

The largest survey for outflows in Perseus to date was un- 
dertaken by Hatchell and coworkers (HFR, Hatchell & Dunham 
2009). They surveyed 83 star-forming cores (compared to our 65) 
for CO / = 3 — » 2 outflows. Our survey is highly complementary 
to the Hatchell et al. survey since we focus on outflow properties 
whereas they concentrate on source classification. Furthermore, we 
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Figure 2. Overview of the new 12 CO J = 3 ->• 2 data towards NGC 1333. 
The colour-scale is integrated intensity, / T£dv from —24 to 40 km s _1 with 
contours at 20, 30, 40, 50, 70, 90, 120, 150, 180Kkms _1 . 



map larger, contiguous areas which encompass a greater amount of 
each individual outflow's emission, thereby increasing the accuracy 
with which we can estimate various outflow parameters. In addi- 
tion, we use our complementary data from 12 CO's isotopologues 
(Paperjj). We correct our outflow mass estimates for the 12 CO op- 
tical depth which we calculate from our 13 CO 7 = 3—^2 data and 
derive estimates of the velocity of the driving source/ambient emis- 
sion from our C 18 data of the same transition. 



3 OUTFLOW DETECTION RATES AND SOURCE 
CLASSIFICATION 

Submillimetre continuum surveys select the earliest stages of star 
formation, when cores are identified as being either starless or 
young protostars. An IR detection of a core or the presence of an 
outflow are thought to be robust indicators of a protostellar source. 
The star-forming cores located by Hatchel l et al.| {2005) with 
SCUBA were subsequently classified from their SEDs (including 
Spitzer data, Hatch elTet al.|2007a) and by the presence/absence of 
outflows (HFR; Hatchell & Dunham 2009 ). A major conclusion of 
HFR is that outflow mapping is as good as Spitzer at identifying 
protostars. 

We use an objective criterion to search for outflows, as ap- 
plied by |HFR| and |Hatchell & Dunham| {2009): are the 12 CO linew- 
ings greater in intensity than 3<7 rms at velocities ±3kms _1 from 
the ambient cloud/driving-source velocity? This definition of de- 
tectable outflowing gas was chosen as it is simple, unbiased, uses 
our ancillary information and allows direct comparison with the 
|Hatchell et al.| surveys. Often, authors rely on more subjective cri- 
teria to identify outflows, such as a definite bipolar structure. When 
we look at measured outflow properties in Section [5] we too will 
only examine outflows with definite extended structure or bipolar- 
ity to ensure we are measuring true outflow characteristics. How- 



ever, the aim of this section is to look in an unbiased manner at the 
classification of SCUBA cores. Our criterion may find a number of 
false positives, due to multiple velocity components along the line 
of sight or other unusual velocity structures in the ambient cloud, 
but we have tried to minimize such contamination by inspecting ev- 
ery spectrum and point out cores with dubious detections in Table 

in 

As the noise on our HARP data varies, we apply this crite- 
rion at ev ery spat ial point individually. We use our existingveloc- 
ity maps {Paper 1 1, computed by fitting a Gaussian functiorjjto the 
C 18 / = 3 — > 2 line profile at every spatial point, as an estimate 
of the ambient cloud velocity, vq. For 11 of the 65 cores surveyed 
(10 cores in NGC 1333 and one in L1455) there was no estimate 
of vo at the same spatial position as the core peak. This is because 
our 12 CO data cover a larger area than our C 18 data in NGC 1333 
and the C 18 line is too weak in one of the cores in LI 45 5. For 
these sources we take vo to be a characteristic velocity for the re- 
gion, listed in Table[T] Therefore we identify blue(red)-shifted out- 
flowing gas at a position if the 12 CO data is greater than 3cr rms 
(where cr rms is calculated for 12 CO at each position) at velocities 
v — (+)3kms _1 . Almost the entire map around NGC 1333 sat- 
isfies this criterion so we increased the offset here to 4kms _1 to 
produce more plausible outflow identifications. 

Table [2] lists the results of the outflow search towards the 
SCUBA cores identified by |Hatchell et al.| ( |2007a] ). Spectra from 
each source are presented in Appendix|B| We record an outflow de- 
tection if all the positions surrounding a source meet the above cri- 
terion. If there is the possibility of confusion with other flows or the 
outflow structure is unrealistic (with small lobes that do not resem- 
ble outflows) we mark it as a marginal detection unless the source is 
particularly well-known. Our fields cover 65 Hatchell et al. (2007a) 
cores: 26 starless, 25 Class and 15 Class I protostars. We detect 
outflows towards 45 (69 per cent) of these (21 firm and 24 marginal 
detections, see Table [3). 

Nearly all of our sources (61 out of 65) were also examined by 
either HFR or Hatchell & Dunham ( 2009). We disagree in our out- 
flow identifications towards 12 cores (18 per cent). Of these twelve 
cores, rsj 7 of our disagreements result from our more conservative 
definition of an outflow, i.e. it has to resemble a standard bow shock 
shape and/or be bipolar for us to record a definite over marginal de- 
tection. Of the remaining five cores (all in NGC 1333), at least three 
have outflowing gas in the vicinity, which we attribute to other driv- 
ing sources. 

All but one of our definite outflow detections are protostellar. 
The exception is HRF85, in the south of NGC 1333 (see Fig.[AT]), 
along the line of blue- shifted gas from IRAS 2 A via HH13. The 
outflow structure is not highly collimated and bipolar, anticipated 
from the source's young age and is therefore perhaps more likely to 
be emanating from IRAS2A. Hatchell & Dunham ( 2009) also note 
the unclear CO outflow structure around this source which makes 
its classification ambiguous. 

Around three quarters of the definite detections are Class 
sources, which is not unexpected considering the greater degree of 
collimation in young outflows (e.g. |Arce & Sargent 2006). Just five 
of the catalogued protostars (all in NGC 1333) do not have any de- 
tection. HRF61 is a prominent non-detection for HFR (and again in 



3 A single Gaussian was fitted as the C ls O J = 3 — >• 2 line has a single 
peak over the majority of the map where it is well detected. In Paper III we 
note only 2 out of 58 spectra, inspected towards the peak of known SCUBA 
cores, have double peaked C ls O/ = 3^2 line profiles. 
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Table 2. Results of the outflow search towards the submillimetre cores in our fields (Hatchell et al. 2007a I. 



Source a 


Region 


Name 


Class b 


Outflow? c 


Confusion? d 


Verdict 6 


Hatchell et al. verdict f 










R 


B 






B 


HARP 


HRF41 


NGC 1 333 


IRAS4A 





y 


y 




Out 


y 


— 


HRF42 


NGC 1 333 


TR A S4R 





y 


y 




Out 


y 


- 


HRF43 


NGC 1 333 


SVS 13 


I 


y 


y 




Out 


y 


— 


HRF44 


NGC 1 333 


TRAS2A 





y 


y 




Out 


y 


- 


HRF45 


NGC 1 333 


SK24 


I 


y 


y 




Out 


y 




HRF46 


NGC 1333 


SK20/21 





y 


y 




Out 


y 


- 


HRF47 


NGC 1 333 


SK31 





y 


y 




Out 


y 


- 


HRF48 


NGC 1 333 


TR A S4C 





y 


n 


y 


M 


y?c41,42 


- 


HRF49 


NGC 1 333 


SK6 


I 


y 


y 




Out 


— 


y 


HRF50 


NGC 1 333 


SK15 


I 


y 


y 


y 


M 


y?c43,51 


- 


HRF51 


NGC 1 333 


SK16 




y 


y 


y 


M 


y?c43,50 


— 


HRF52 


NGC 1 333 


SK14 





y 


y 


y 


M 


y?c43,50 


— 


HRF53 


NGC 1 333 


SK26 


s 


y 


y 


y 


M 


y?c45 


— 


HRF54 


NGC 1 333 


SK28 


I 


y 


n 




M 


— 


y?c56 


HRF55 


NGC 1 333 







n 


n 




No 


- 


y?c45 


HRF56 


NGC 1 333 


SK29 


I 


y 


y 


y 


M 


- 


y 


HRF57 


NGC 1 333 


SK33 




n 


n 




No 


— 


n 


HRF59 


NGC 1 333 






y 


y? 


y 


M 


- 


y?c43 


HRF60 


NGC 1 333 






n 


n 




No 




y?c45 


HRF61 


NGC 1 333 







n 


n 




No 


11 


n 


HRF62 


NGC 1 333 


SK18 





n 


y 


y 


M 


y?c46 


— 


HRF63 


NGC 1 333 


SK32 


I 


y 


y? 




M 




y? 


HRF64 


NGC 1 333 






n 


n 




No 


n 


n 


HRF65 


NGC 1333 


SKI 





y 


y 


n 


Out 




y 


HRF66 


NGC 1333 


SK30 


S 


y 


y 


y 


M 


- 


y?c56 


HRF67 


NGC 1333 




I 


y 


y? 


n 


M 


- 


y 


HRF68 


NGC 1333 







y 


y 


y 


M 


y?c45 




HRF69 


NGC 1333 




I 


n 


y 


n 


M 


- 


y?c44 


HRF70 


NGC 1333 


SK22 





n 


y? 


n 


No 


- 


y?c46 


HRF71 


NGC 1333 







y 


n 


n 


M 


- 


y? 


HRF72 


NGC 1333 




S 


n 


n 


n 


No 


- 


y?c43 


HRF74 


NGC 1333 




I 


n 


n 


n 


No 


- 


y 


HRF75 


NGC 1333 







n 


n 


n 


No 


- 


y 


HRF85 


NGC 1333 




S 


y 


y 


n 


Out 


- 


y? 


Bolo26 


NGC 1333 




S 


n 


n 


n 


No 


- 


y?c65 


Bolo44 


NGC 1333 




S 


y 


y 


y 


M 


- 


- 


HRF12 


IC348 


HH21 1 





y 


y 




Out 


y 


y 


HRF13 


IC348 


TC34R-MMS 





y 


y 




Out 


y 


- 


HRF14 


IC348 




I 


y 


y 




Out 


y 


- 


HRF15 


IC348 


TC34R-SMM3 





y 


y 




Out 


y 


y 


ni\r 1 o 






c 



y? 


n 


n 


1V1 


n 


n 


HRF17 


IC348 




S 


y? 


y? 


y 


M 


y?cl3? 




HRF18 


IC348 




S 


y? 


n 


n 


No 


n 


n 


HRF19 


IC348 




S 


n 


n 


n 


No 


n 


n 


HRF20 


IC348 




S 


n 


n 


n 


No 


n 


n 


HRF21 


IC348 




S 


y? 


n 


n 


No 


n 


n 


HRF23 


IC348 




S 


n 


n 


n 


No 


n 


n 


HRF24 


IC348 




S 


n 


n 


n 


No 


n 


n 


HRF25 


IC348 




S 


n 


n 


n 


No 


n 


n 


HRF26 


IC348 




S 


n 


n 


n 


No 


n 




HRF101 


IC348 




I 


n 


y 


n 


M 




y 


Bololll 


IC348 




S 


n 


n 


n 


No 






Bololl3 


IC348 




S 


n 


n 


n 


No 







a Hatchell et al. ( 2007a I source number. 

b Source classification from |Hatchell et al.|{2007a) : S=starless, 0=Class and I=Class I. 

c Is red (R) or blue (B) outflowing gas detected according to the outflow criterion (see text)? Sources are marked '?' if only a limited region has outflowing 
gas. 

d Is the outflow confused with others? 

e Out=definite detection, M=marginal and No=no detection. 

Outflow verdict of the Hatchell et al. surveys using either receiver B |HFR| or HARP jHatchell & Dunham 2009) on the JCMT. Question marks denote 
potentially confused sources, with any numbers recording the sources of such confusion. 
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Table 2 - continued 



Source 


Region 


Name 


Class 


Outflow? 


Confusion? 


Verdict 


Hatchell et al. verdict 










R 


B 






B 


HARP 


HRF27 


L1448 


L1448NW 





y 


y 


y 


Out 


y 




HRF28 


L1448 


L1448N:A/B 





y 


y 


y 


Out 


y 




HRF29 


L1448 


L1448C 





y 


y 


n 


Out 


y 




HRF30 


L1448 


IRS2 





y 


y 


n 


Out 


y 




HRF31 


L1448 







y 


y 


y 


Out 


y?c30? 




HRF32 


L1448 




S 


y 


n 


y 


M 




n?c27 


Bololl 


L1448 




S 


y 


n 


y 


M 






HRF35 


L1455 


RN015-FIR 


I 


y 


y 


n 


Out 


y 




HRF36 


L1455 







y 


y 


n 


Out 


y 




HRF37 


L1455 




I 


y 


y 


y 


M 


y? 




HRF39 


L1455 




I 


y 


y 


y 


M 


y? 




HRF40 


L1455 




s 


y 


n 


y 


M 


y?c35,36 





HARP data, Hatchell & Dunham 2009), which they explain might 
be caused by radiation and winds from nearby luminous stars blow- 
ing away its surrounding gas. The other non-detections in HRF55, 
HRF70, HRF74 and HRF75 might be explained similarly or by 
having extremely weak and/or small flows. Outflows close to the 
plane of the sky will also be hard to detect. Furthermore, orienta- 
tion and a non-uniform environment possibly explain the protostars 
with only one prominent outflow lobe. 

Of the marginal outflow detections towards starless cores, all 
but one (HRF85) are potentially confused with strong outflows 
from neighbouring protostars. Therefore we exclude them from the 
following analysis. Our outflow detection rates are comparable to 
HFR s; the extra sensitivity of this survey (over HFR ; compare their 
sensitivity of ~ 0.3 K in lkms" 1 toTable[j} has not overturned any 
of the starless classifications of Hatchell et al. ( 2007a) as found by 
|Hatchell & Dunham| ( |2009t as well. Indeed, we are probably lim- 
ited more by resolution - many of the single SCUBA clumps are 
well-known binaries with multiple outflows not resolved with the 
JCMT. Additionally, confusion (due to resolution and crowding) is 
a severe hindrance in classifying sources in this way. The original 
continuum survey is incomplete to low-mass protostars: Hatchell 
|et al.| p007a| ) calculate the 3(7 completeness limit is 0.3 M for 
a 10 K protostar. This estimate of the mass limit is uncertain by 
up to a factor of ~ 10, given the uncertainties in the properties of 
the dust. By searching for outflows only towards identified SCUBA 
cores we miss outflows driven by low-mass protostars, which were 
undetected in the continuum survey. 



4 ESTIMATION OF OUTFLOW PARAMETERS 

In this section, we describe the methods we will subsequently em- 
ploy to compute the physical parameters of extended outflows, 
identified around the protostars in our survey fields in Sectionp] 
These properties for individual outflows are detailed in Tables |B1| 
and |B 2 1 with the totals listed for all the outflows in each region in 
TableS 



4.1 Mass 



We follow a standard method (e.g. Goldsmith et al. 1984) and cor- 
rect for the optical depth of 12 CO using our complementary 13 CO 
/ = 3 — > 2 data. At moderate outflow velocities (close to that of 



the ambient cloud emission), 12 CO is optically thick but any de- 
tectable 13 CO is optically thin. Assuming 12 CO is optically thick, 
13 CO optically thin, local thermodynamic equilibrium (LTE), iden- 
tical beam-filling factors and excitation temperatures for both iso- 
topologues, the following relation holds (e.g. Hatc hell et al.|i 999 ): 



r;(i 3 co) 



V12V 

Vl3/ 



[ 12 CO] 1 - 



-T12 



[ 13 CO] T12 



(1) 



where V 12 = 345.796GHz and V 13 = 330.558 GHz are the tran- 
sition frequencies of the 12 CO and 13 CO lines respectively, T12 
is the optical depth of the 12 CO gas and t he abundance ratio, 
[ 12 CO]/[ 13 CO] = 62 ( |Langer & Penzias|l993| ). z u can be evalu- 
ated numerically and the correction, — exp(— T12)), applied 
to the 12 CO mass. Cabrit & Bertout ( 1990 ) showed that this method 



was sufficient to correct the mass for optical depth effects in typ- 
ical outflows. We apply the correction to all the velocity channels 
where there is detectable 13 CO emission above the noise (i.e. where 
it is > c rms ). The 13 CO data do not cover the entire area observed 
in 12 CO for NGC 1333, so some of the outflow masses are not ad- 
justed. A majority of outflows are corrected for the optical depth 
of 12 CO (28 out of the 32 outflows whose parameters we report in 
Tables [bT] and |B2[ ). The masses are increased by a factor ranging 
from 1.8 to 14.3 with a median of 3.8. 

Without an opacity correction, the column density of CO, 
Af(CO), can be calculated from the corresponding / = 3 — » 2 emis- 
sion assuming LTE thus: 

^ggn = 4.72xio" , wk,„ ^/w^y (2) 



V cm- 2 J A " exp(-33.7K/r ex ) V Kkms" 1 

where T ex is the excitation temperature, which following HFRlwe 
take as 50 K. For an outflow in Perseus (250 pc away) containing 
Npix = 100 pixels, each of area A p i x and average integrated antenna 
temperature, (/T^(v)dv) = 1 Kkms -1 , its mass is: 

1 



(£)-»«^(£0" m 



^pix 



100 



9 sq. arc sec 
lKkms- 1 



(3) 



using the value of the abundance of CO relative to H2, Xqo = 10 -4 , 
from Frerking, Langer & Wilson (1982). In our analysis we inte- 
grate the emission to find the mass over the total observable veloc- 
ity extent of the outflow, i.e. from 2 km s _1 away from the ambient 
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Table 3. Number of outflow detections towards SCUBA cores. 



Core type 


Outflow 


Marginal 


Non-detection 


Total 


Starless 


1 


10 


15 


26 


Class 


15 


5 


4 


24 


Class I 


5 


9 


1 


15 


Total 


21 


24 


20 


65 



cloud velocity to the velocity offset, above which the outflow emis- 
sion falls to below the noise (v max ; see also Sections |4~2] and |43] >. 

In addition to any uncertainties associated with the above 
equations, masses derived from CO observations may not repre- 
sent the true mass of outflowing gas for a number of reasons: (i) 
Strong shocks can dissociate molecular gas into atomic material 
(e.g. Downes & Ray 1999 ) so we may miss a significant fraction of 
the total mass. Fortunately, outflow mass estimates are dominated 
by low-velocity material which is unlikely to have passed through 
strong shocks; (ii) Some of the low- velocity outflow emission may 
be missed as it is swamped by signal from the ambient cloud mate- 
rial. If we consider material moving at a modest offset (< 2 km s _1 ) 
from the ambient cloud velocity to belong to the outflow, |Downes| 
& Cabrit (2007 ) estimate that the mass underestimate could be as 
large as a factor of 2-3; and (iii) In reality outflows will not be at 
a single constant temperature. Downes & Cabrit (2007 ) show that 
assuming a temperature of 10 K for mass estimates from the CO 
/ = 2 — >> 1 transition introduces only a small error into the outflow 
momentum and none into the mass. We expect similar results with 
our higher assumed temperature for the / = 3 — >> 2 line. 



4.2 Momentum and energy 

We define the observed outflow momentum along the jet axis: 
J m(v)|v — vq| dv, 



Pout 



(4) 



where vo is the outflow centre velocity, approximately the velocity 
of the driving source, which we derive from our C 18 07 = 3 — >>2 
observations, and m(v) is the observed 12 CO mass corrected for 
the 12 CO optical depth. Inclination effects can greatly reduce the 
amount of momentum observed (by 1 / cos i for uni-directional jets, 
where i is the angle of inclination to the line of sight). We do not 
make any inclination corrections until Section [578] as the inclina- 
tion angles are mostly unknown. The largest contribution to the 
momentum comes from the highest- velocity material, that most af- 
fected by dissociation. Therefore, any outflow momentum derived 
from CO observations further underestimates the true momentum, 
with the estimate getting poorer as the density contrast between the 
jet and ambient material increases (Downes & Cabri t|2007] >. 
An outflow's kinetic energy is similarly: 



\f m[y)\\ 



-v | 2 dv. 



(5) 



This is even more affected by inclination (1 / cos 2 i correction) and 
dissociation. Downes & Cabrit ( 2007} also point out that the con- 
stant temperature assumption will underestimate the kinetic energy 
for an outflow with high density contrast, which can only be recti- 
fied by assuming a higher excitation temperature for the fast mate- 
rial over the slow. 

We integrate over the full observable velocity-extent of the 
outflow, as above for the mass calculation (see Section [4T| . This 
yields the best estimates of individual outflow properties and allows 



us to examine more accurately the total impact on their natal envi- 
ronment. As the maximum velocity observed depends on the noise 
level, a larger fraction of less massive outflows will be missed. An 
alternative method, employed by HFR reduces this bias by remov- 
ing the high- velocity contribution to all the outflows through only 
integrating between set velocity offsets. However, depending on the 
velocity limits chosen, this method might severely under-estimate 
the outflow momentum and energy as the largest contribution to 
these comes from the highest outflow velocities. 



4.3 Dynamical time and driving force 

The outflow momentum flux or driving force/thrust, Fqq, is a cru- 
cial input into outflow models and the parameter used to infer that 
outflows decrease in strength with age (Bontem ps et al.| [l996). 
Methods to estimate Fqq typically fall into two categories: 'global' 
or 'local'. Local measures examine the outflow in the vicinity of 
the protostar, which is the only option if the maps do not cover all 
of the outflow. We use a global measure, looking at the entire out- 
flow, as our datasets are sensitive and extend over a large portion 
of the observable outflows. Downes & Cabrit (2007 ) evaluated a 
number of methods to find Fqq using a jet-driven outflow model 
(most suitable for Class protostars which are the majority of our 
driving sources) and we follow the route they find to be most ac- 
curate: the so-called 'v max ' method <Lada|19 85; Cabr it & Bertout| 
|19921|Beuther et al.|2 002). We define the dynamical time, T d , as the 
time for the bow shock travelling at the maximum velocity in the 
flow, v max , to travel the projected lobe length, Li be : 



Aobe 
v max 



(6) 



We estimate v max as the maximum velocity observed across the 
outflow lobe (the manner in which we assign spectra to individ- 
ual outflows is described in Section \5A\ i.e. the maximum velocity 
above which the outflow emission falls below the noise, Oj-mS' 

The 

force can then be estimated: 



^co 



Pout 



(7) 



The main uncertainty comes from t<j. If the expansion rate of the 
outflow varies, any estimate of is inaccurate without detailed 
knowledge of the mass-loss history of the protostar. Parker, TacT] 
|man & Scott| (l991) suggest may underestimate true outflow 
ages by up to an order of magnitude. This may well be true for 
evolved sources where the flow has left the parent cloud but not 
young Class sources, where Downes & Cabrit (2007 ) suggest the 
outflow age is overestimated by the dynamical time-scale. Due to 
inclination effects, even the maximum outflow velocity does not 
necessarily reflect the true bow advance speed. Such inclination 
effects are worst for in outflows close to the plane of the sky, 
where lobes appear their longest and v max is its smallest. However, 
we cannot account for this without accurate inclination angles (cal- 
culated from e.g. proper motion studies) for a large number of the 
outflows. 



5 OUTFLOW AND PROTOSTELLAR EVOLUTION 
5.1 Parameters and outflow allocation 

In Section [3] we extracted the CO J = 3 — > 2 spectrum at the peak 
position of every SCUBA core jHatchell et al.|2005| [2007a] in our 
fields and looked for an outflow. This method produced similar 
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starless/protostellar core classifications as those based on Spitzer 
SEDs (see also |HFR| |Hatchell & Dunham||2009| >. For the rest of 
this paper we trust that the SED classifications are accurate and 
now look for extended outflow structure around the 34 identified 
protostars in our fields ( Hatchell et al. 2007a ). We decide if a spa- 
tial position has blue- or red- shifted outflows based on the objec- 
tive criterion detailed in Section [3] is the intensity in the linew- 
ings > 3(7 rms ±3 km s _1 from the the driving- source velocity? Con- 
tiguous regions adjacent/overlapping aprotostar, whose constituent 
spatial positions satisfy this criterion are the ones included in an 
outflow's blue or red lobe. For reasonably isolated protostars this 
easily defines the outflow structure, however, in more crowded re- 
gions (particularly NGC 1333) the outflows from different sources 
overlap and are confused (see Table [2). For confused sources we 
allocate the emission, selected with the outflow criterion, to differ- 
ent source lobes more subjectively. We were guided by integrated 
intensity contours of the outflows (see Figs. |Al| to |A9| > and con- 
straints from any corresponding red or blue lobe on the opposite 
side of the source, which are often not confused (e.g. we expect 
these to be in a straight line). 

For the outflow lobes identified, we estimate various physi- 
cal parameters using the prescriptions in the previous Section ([4]). 
The physical properties of the outflows are listed for each source 
in Tables |Bl| and |B2| of Appendix [B] In addition, Table [^summa- 
rizes the different trends in parameters for the outflow population. 
In general the parameters we derive are comparable to other single- 
dish surveys (see pavis et al.|l 997 ||KS00l|Taf alia et al.|2006) but 
slightly larger. Small discrepancies are unsurprising as the alloca- 
tion of emission to particular outflows tends to be rather subjective, 
particularly in confused regions. Our parameters are consistently 
slightly higher than previous estimates as we: (i) apply a correction 
for the 12 CO optical depth; and (ii) have increased sensitivity and 
areal coverage in our maps. However, the dynamical time-scales 
we estimate are often less than half of estimates using the intensity- 
weighted velocity. The v max method should always produce com- 
paratively smaller time-scales. This tendency is reinforced in our 
estimates as we are able to follow the outflow emission to higher 
velocities because our observations are more sensitive. is par- 
ticularly inaccurate for older Class I outflows, which are expected 
to be longer and wider. Many of their terminal bow shocks may 
have left the gaseous cloud, so cannot be the outflow age. Addi- 
tionally, longer, older flows have had greater opportunity to interact 
with other outflows and obstacles, so they are more confused. As 
an illustration, SVS13, a Class I source, drives a fast, confused out- 
flow. Using v max , we find = 4380 yr, younger than many of the 
Class flows, whereas [KS00| estimate T dyn = 30510 yr. SVS13 is 
an extreme case but demonstrates the inherent uncertainty in the 
derivation of and therefore Fco 



5.2 Length 

A simple expectation is that older flows are longer and those in 
crowded regions appear shorter as they intersect with others. There 
is little difference in the mean length between blue- and red- shifted 
lobes, being (120±2oQand (130 ±20) arcsec respectively. For 
any particular outflow, the lobes do not extend symmetrically ei- 
ther side of the driving source. The lobe length ratio (the larger 



4 All the uncertainties quoted are errors on the mean, cr m formed from the 
standard deviation of the values a and their number, N, via cr m = a/y/N. 



over smaller lobe length) is on average 2.0 ±0.6 with no signifi- 
cant differences between sources of different types or region. 

Outflows appear to get shorter during source evolution, aver- 
age lobe lengths are (140 ±20) and (88 ± 11) arcsec for Class 
and I sources respectively. There is a lot of scatter but it is signifi- 
cant that the longest outflows are all Class (see Table |B1| . How- 
ever, the apparent length of an outflow is not necessarily related to 
its age. This obviously has major implications for other parameter 
estimates such as and Fqo- Flows from Class sources are typ- 
ically found to be stronger and more collimated than those from 
Class Is ( |Bontemps et al.|1996[|Arce & Sargent|2006| >. Indeed, on 
further evolution into the Class II stage, outflows seem to lose all 
definite structure (Arce & Sargent 2006 ). We already noted that it 
is more difficult to identify Class I outflows: three quarters of the 
definite detections were Class 0. Presumably the same is true for 
the length of a Class I flow; if material gets more spread out with 
time, at lower column density, then we will not necessarily detect 
the weaker parts which may shorten its apparent length. In addi- 
tion, it is more likely that older Class I outflows will have broken 
free from the parent cloud, so we would underestimate their full 
extent. 

The different regions also seem to have an effect on the 
observed length, although there are small numbers of outflows 
everywhere apart from NGC 1333. The longest outflows are in 
L1448, (240 ±40) arcsec, with NGC 1333 and L1455's under half 
their length, at (110 ±20) and (90 ±30) arcsec respectively while 
those in IC348 are the shortest: (50 ± 9) arcsec. L1448 has ex- 
tremely long outflows, extending perpendicular to the molecular 
ridge. |Bally et al.| ( |2008) ) suggest that the outflows from L1448 and 
LI 455 are huge, stretching towards each other and possibly inter- 
acting. All the outflow sources in L1448 are Class 0, so their long 
lengths greatly affect the overall region average. The outflows in 
NGC 1333, which could presumably stretch that far if undisturbed, 
are probably appearing shorter as they overlap. 



5.3 Maximum observed outflow velocity 

The trends in v max affect all the derived outflow parameters. Class 
sources have higher- velocity outflows; on average v max is (17.9 ± 
1 .4) km s -1 for Class protostars and (12.2± 1.6) km s _1 for Class 
Is. This supports the idea that Class outflows are stronger. Inter- 
estingly, the trends in lobe length and v max will counter one another 
in Td- The longer length of Class outflows increases T^, while 
their higher velocity decreases T^. Overall, this may explain a lack 
of differentiation between for different outflow stages. 



5.4 Mass, momentum and kinetic energy 

The derived outflow masses range from 0.002 M for that from 
HRF71, a Class I source, to 0.4 M for the vast Class outflow 
from NGC 1333-IRAS2A. In general, Class outflows are more 
massive than Class Is (see Fig.[3j, hardly a surprise since we found 
that Class outflows appear longer (see Section \5.2\ and span 
larger velocity ranges (Section [53] >. On average, the masses are 
(Moat) = (0.09 ± 0.02) M for Class and (0.06 ± 0.03) M for 
Class I outflows. The distributions are not significantly different, a 
Kolmogorov-Smirnov (K-S) test is inconclusive with a 9 per cent 
probability that they are drawn from the same underlying popula- 
tion. Across the different regions, the mean mass is not very reveal- 
ing, with most populations severely affected by the small number 
of sources. 
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Table 4. Summary of the trends in physical outflow parameters as defined in the text. The outflows are grouped by region or by the classification (0 or I) of 
their associated driving source. All the errors listed are errors on the mean not sample deviations. 



Population 




<5 L 


Vmax 


c v 


M out 




Pout 


G P 


^out 


o E 






^co 






(arcsec) 


(kms~ 


- 1 ) 


(M ) 


(M© 


kms" 1 ) 


(10 36 J) 


(yr) 




(10" 5 


M© kms -1 yr -1 ) 


NGC1333 


110 


30 


15.3 


1.4 


0.09 


0.03 


0.5 


0.2 


9 


4 


9800 


1500 


7.6 


3.7 


IC348 


50 


9 


12 


3 


0.009 


0.005 


0.05 


0.02 


0.7 


0.4 


5000 


500 


0.9 


0.4 


L1448 


240 


40 


24 


3 


0.15 


0.03 


1.4 


0.8 


36 


17 


14000 


5000 


16.7 


11.4 


L1455 


90 


30 


11.5 


0.6 


0.032 


0.006 


0.15 


0.03 


1.7 


0.3 


10000 


3000 


1.9 


0.7 


Class 


140 


20 


17.9 


1.4 


0.09 


0.02 


0.7 


0.2 


14 


5 


9800 


1600 


7.8 


3.0 


Class I 


88 


11 


12.2 


1.6 


0.06 


0.03 


0.3 


0.2 


7 


6 


10000 


2000 


6.8 


5.8 


All 


124 


14 


15.9 


1.1 


0.08 


0.02 


0.54 


0.17 


11 


4 


9800 


1200 


7.4 


2.8 



a Lobe length. 



The outflow momenta range from 0.007 M kms 1 for 
HRF15, the small bipolar outflow from IC348-SMM3, to 
4.4 M kms" 1 for HRF30, one of the energetic outflows in LI 448. 
These two outflows also have the lowest and highest energies: 
5 x 10 34 J and 1 x 10 38 J respectively. We plot histograms of the 
momenta and energies in Fig. [3] The pattern for both parameters is 
nearly identical. High-velocity outflow material makes an increased 
contribution to the outflow energy over momentum (Downes & 
Cabrit 2007), assuming a power-law velocity dependence of the 



mass, m(v) « v ^, where 1.5 ^ y ^ 2.0 (e.g. |Smith, Suttner & 



|Yorke|1997[|Downes & Rayp 999 ). A small difference in the pat 
tern of values between momentum and energy implies there is little 
high- velocity material in the flow i.e. these are unusually steep out- 
flows with high values of y or alternatively the outflows span only 
a limited velocity range. Given that the maximum detected outflow 
velocity is typically 15 kms -1 , the former is probably favoured. 
The largest Class I momentum and energy bins in the histograms 
both only contain SVS13, with the second largest bins empty (the 
largest mass bin also just SVS13), suggesting SVS13 is a highly- 
powerful Class I anomaly. However, it is a particularly compli- 
cated source and is unlikely to harbour a single Class I protostar 
(see Section [A 1.4) . HFR suggest its outflow may be driven by a 
Class protobinary companion. With SVS13 included, the outflow 
momentum is marginally higher for Class than Class I sources: 
(Pout) = (0.1 ± 0.2) and (0.3 ± 0.2) M© km s _1 respectively. If we 
leave out SVS13, the Class I momentum lowers considerably to 
(0.10 ± 0.03) M km s _1 , making the Class I population have sig- 
nificantly less momentum than the Class 0. A K-S test comparing 
the two samples (without SVS13) yields a lper cent probability 
that they are drawn from the same population. The energy has a 
similar dependence on SVS13 but is more extreme. Without SVS13 
the Class population has on average over ten times more kinetic 
energy than the Class I: (E out ) = (1.4 ±0.5) x 10 37 compared to 
(1.0 ± 0.3) x 10 36 J respectively. There is then an even smaller K-S 
probability (0.4 per cent) that they are drawn from the same popu- 
lation. 



5.5 Dynamical time-scale 

The dynamical time- scale, does not show the same differentia- 
tion between the source populations (see Fig. |4j, with averages of 
(9800 ± 1600) and (10000 ± 2000) yr for Class and I outflows 
respectively. Indeed, a K-S test finds an 89 per cent probability that 
they are drawn from the same population. This is not a surprise: 



Class outflows are longer and faster, effects that act against each 
other to produce similar average to Class Is. 

In Fig. [4] we compare to the age of the driving source 
The bolometric temperature, T^oi, is the temperature of a black- 
body with the same mean frequency as the protostar' s SED'7], i 
should increase as a protostar clears its circumstellar dust and pro- 
gresses towards the main sequence ( Mye rs & Ladd|1993]>. It is thus 
a measure of the source age, £yso Hatchell et al. ( 2007a ) use it to 
differentiate between Class and I protostars, with the former be- 
low and latter above a division at T^x = 70 K ( |Chen et aT7f l995 ). 
|Ladd, Fuller & Deane] {1998 ) find a best-fitting relation between 
*yso and r bol : 



log(Wyr) = [2.41og(r bol /K) -0.9] ±0.6, 



(8) 



which we use to compute tyso f° r eacn driving source, taking T^ \ 
from Hatchell et al. (2007a). Using this relation the break between 
classifications occurs at tyso — 3375 yr. There is no correlation be- 
tween Td and tyso for either population or overall (see Fig. [4). 

The dynamical time is not the age of an outflow (Padman, 
Bence & Richer 1997}. Statistical outflow lifetimes deduced by 
|Parker et al.| ( |199lt and |Fukui et al.| ( |1993| ) are at least an order of 
magnitude larger. depends on the time history of the jet velocity 
and the density of the region between the driving source and ter- 
minal bow shock. This assumes we can identify the terminal bow 
shock, i.e. it is not confused or lies beyond the map boundary. A 
simple calculation for the oldest sources in this sample with ages 
of around 10 6 yr, implies an outflow length of ~150pc for a jet ve- 
locity of 150 km s _1 - far larger than the area we mapped. Even 
for Class sources, whose shorter, brighter and more collimated 
outflows are easier to accurately trace, there is little correlation be- 
tween Td and age. 

In summary, is a poor estimate of the outflow age, with 
little differentiation between Class and I outflows. However, the 
problem for parameters derived from may actually not be as se- 
rious. Presumably we may underestimate the outflow age from 
if we do not detect the terminal bow shock or do not trace the out- 
flow's extent accurately because of confusion or low column den- 
sities of material. In such a case we also underestimate the total 
mass, momentum and energy. To some extent these underestimates 



5 Implicit in all the analysis that follows, is that the outflow and protostellar 
ages are the same, i.e. outflows endure for the entire embedded phase, which 
is a fairly sound deduction from the high outflow detection rates towards 
protostars. 
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Figure 3. Histograms of outflow mass (top), momentum (middle) and en- 
ergy (bottom). The dashed light blue line is for the Class outflows, whilst 
the solid purple one is for Class Is. The largest mass, momentum and en- 
ergy bins for the Class Is only contain the outflow from SVS13. The bars 
represent Poisson errors on the bins. 



cancel each other out. Other methods that characterize an outflow 
by its emission close to the driving source may in fact be a fairer 
way of estimating parameters (e.g. Bontemps et al. 1996, HFR) 
where there is less confusion and possibility of missing emission, 
although we have already noted that some of these methods may 
underestimate the flux in an entire outflow. 



5.6 Momentum flux 



Bontemps et al.|(T996| ) measured the momentum flux, Fqq, for a 
sample of 9 Class and 36 Class I outflows. They found the mo- 
mentum flux declined with age, implying a similar decrease in the 
mass accretion rate onto the central object (M acc ) with this rate di- 
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Figure 4. Upper panel: Histogram of outflow Td, again the bars repre- 
sent Poisson errors on the bins. Lower panel: Td, versus nominal driving 
source age, ?yso, derived from the bolometric temperatures in Hatchell et al. 
(2007a) . 



rectly proportional to the source envelope mass (M env ). We look for 
a similar trend in Fig. [5] We estimate the flux in two ways using dif- 
ferent time- scales: (i) Fqo uses the conventional (Fqo = Pout/ ^d) 
and (ii) F£ takes the age of the outflow to be the same as the proto- 
star, ^yso (^co = Pout Ayso)- The same trends are apparent for both 
estimates but F£ has greater differentiation between the types. 

Class sources have higher Fqo than Class Is. Most of the 
high Fco sources are Class apart from SVS13 which, as al- 
ready mentioned, has an anomalously large momentum for a Class 
I source. The sample averages bear this out if SVS13 is excluded 
from the Class I population: (F co ) = (0.8 ± 0.3) x 10" 4 for Class 
outflows and (1.1 ±0.3) x 10~ 5 M krns" 1 yr" 1 for Class Is 
((0.7 ±0.6) x 10" 4 with SVS13 included). F* Q accentuates the 
same trends; the Class I average is over an order of magnitude 
smaller than the Class 0: (6 ±2) x 10~ 6 ((1.3 ±0.8) x 10~ 5 with 
SVS 13 included) compared to (4.3 ±1.5) x 10 _4 M© kms" 1 yr _1 . 
For both Fco and F£ there is an extremely small K-S probabil- 
ity that the two protostellar classifications are drawn from the same 
population: < 0. 1 per cent. Thus, the outflow momentum flux de- 
clines from the Class to I stage as found by Bontem p*s et al.| In- 
triguingly, HFR do not find this relation holds. Their results are 
influenced by a small number (3-4) of high-flux Class I and low- 
flux Class outflows. The former they explain are potential bi- 
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Figure 5. Histograms of the outflow momentum flux. The momentum flux 
is derived using a time-scale of Td C^co, upper panel) or the age from 7t, i 
(Fq , lower panel). The bars represent Poisson errors on the bins. 



nary sources, containing a more energetic Class companion. Af- 
ter SVS13, HFR find the two most powerful Class I sources are 
HRF45 and HRF50 both in NGC 1333. We also measure a large 
Fqo for HRF45 (the second largest), although it is considerably 
smaller than for SVS13 and close to the other Class I protostars. 
On the other hand, we do not evaluate Fco for an outflow from 
HRF50 as it is so confused with the one from SVS13 (see Section 
[AL4}. 

Low-flux Class sources are harder to explain; perhaps some 
of the youngest Class flows are just powering up and have not 
reached their maximum outflow force (e.g. [Smith 2000). Further- 
more, the source classification based on SEDs will not be entirely 
infallible. Inclination effects for protostars with a disc/bipolar cav- 
ity can cause large deviations from typical SEDs for their class (e.g. 
|Whitney et al.|2003] > causing them to be mis-classified. 

The correlations of the momentum flux (both Fqo and Fq ) 
with various properties are shown in Fig. [6] Since the range of pa- 
rameters is small and the fluxes are not corrected for inclination, 
there may be considerable uncertainty. The first correlation is with 
T^oi and if the outflow force does decrease with age, we expect there 
to be a decrease in Fqo with T\, \, which is tentatively apparent. Of 
course with F£ there is a much firmer trend, expected as F£ by 
definition is inversely related to 7b i- Hogerheijde et al. (1998) also 
note that an outflow's age is not a good predictor of the force (the 
conventional Fqo). 

The second correlation, with bolometric source luminosity, 
Lb i, is well known (e.g. |Cabrit & Bertout| |l992). We take L^ \ 
from Hatc hell et al.| ( |2007a| ), who assumed a different distance to 
Perseus. Correcting for the distance would uniformly reduce L^oi 



by 40 per cent, leaving the shape of the correlation unaffected. 
In these data there seems to be two distinct correlations for the 
different evolutionary stages in both Fco and F£ with slopes 
very similar to unity. Both types have forces many times larger 
than the momentum supplied in radiation from the driving source 
C^rad = Lfcoi/c). The correlations are approximately consistent with 
Class I sources having Fco = 100F ra d with the flux of Class 
sources ~10 times higher again. For Fqo> maximum-likelihood fits 
follow: 



log^co/Mokms-^r- 1 ) = (0.5 ±0.2) log(L bol /L ) 
-(4.8±0.2) Class sources, 



log(F C o/M kms yr 



) = (O.7±O.4)log(L bol /L ) 
-(5.6±0.2) Class I sources. 



(9) 



(10) 



Class objects drive systematically more powerful outflows than 
Class Is for the same luminosity of object, implying they are more 
efficient at driving their outflows. This was also suggested by |Bon-| 
temps et al. ( 1996), although they had fewer Class sources which 
did not allow a good constraint on the correlation. 

The final correlation with M env is the most convincing. Again 
we use M env from Hatch ell et al.| ( |2007a| ), who derive them from 
SCUBA fluxes assuming a distance of 320 pc to Perseus. Like the 
Lboi correlation, to correct for the distance requires a systematic 
reduction in M env by ~40 per cent. Maximum-likelihood fitting to 
the entire population yields a near direct-proportionality relation: 



log(F co /M kms yr 



(O.85±O.16)log(M env /M ) 

-(5.4 ±0.2) (11) 



This correlation again points to a decline in outflow force with 
source age, if the envelope mass is inversely related to the outflow 
age, i.e. Class sources have large masses and Class Is small. 

The common gradient for all sources on a Fco-L^ \ diagram 
led many authors to conclude that there is a common driving mech- 
anism for all outflows and their energetics are dominated by the 
luminosity of the central source (e.g. Lada 1985, Levreault 1988). 
However, as Bontemps et al. ( 1996} showed (and we have seen), 
there are considerable differences in the energetics of outflows at 
different evolutionary stages. To emphasize such trends, Bontemps 
|et al.| tried to remove the luminosity dependence from the Fco-^env 
relation (their fig. 7). We produce an equivalent plot in Fig. [7] fol- 
lowing the appropriate prescription from HFR] Unlike [HFR there 
is a clear evolutionary trend with Fco-c/^bol declining markedly 
from Class to I sources. Thus, we confirm the findings of |Bon-| 
temp s et aT] ( |1996|), and demo nstrate the increased efficiency in 
Class outflows. |Bontemps et al.| parameterize the momentum flux: 



Fco = /er 



M aC( 



-v w xM a( 



(12) 



where / ent is the entrainment efficiency (< 1), i.e. the fraction 
of matter from the central jet/wind fed into the the outflow, M w 
the mass-loss rate in the jet/wind and v w the jet/wind velocity. 
In such a parametrization they state the only likely cause of a 
large decline in Fqo is a reduction in M acc as the other parame- 
ters do not vary enough in plausible jet and X-wind models. How- 
ever, it seems highly unlikely that the entrainment efficiency will 
be the same for Class and I outflows given the differences in 
their surrounding environments and degrees of collimation. Nev- 
ertheless, if we adopt similar values for the parameters to IBon-| 
temps et al. / ent = 1, M w /M acc = 0.1 and v w = 150kms , the 
reduction in Fqo from Class to I sources, (0.8 ±0.3) x 10~ 4 to 
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Figure 6. Correlations of Fco (upper panels) and F£ (lower panels, see text) with r bo i (left), L bo i (middle) and M env (right) taken from Hatchell et al. ( 2007a). 
Sources are marked according to their different Hatchell et al. ( 2007a) classifications. Dashed lines represent the force available in stellar photons, Fco =L\y \/c 



(for a single scattering). Solid lines mark maximum-likelihood fits to the overall population (unless specified): (a)7 7 co 06 T^ 



4 , (b)Fco-A* 



(blue, 



Class 0) and F co ~ L bo i a7 (purple, Class I), (c)F C0 <* M cm °- 9 , (d)F£ °< Z^i" , (e)F£ oc L bo i°- 51 (blue, Class 0) and F£ ~ L bo i°- 8 (purple, Class I) and 
(f)F c * ocM env L7 . 



(I.l±0.3)xl0 -5 Mekms -1 yr 1 on average, implies a reduction 
in the accretion rate from M a cc = 5 x 10 -6 to 7 x 10 -7 M yr -1 , 
somewhat smaller than the reduction inferred by Bontemps et al. 



5.7 Destruction time-scale 

Outflows act as one means of removing matter from a star-forming 
core. The time that it would take an outflow to eject the remain- 
ing mass of its parent core is a limit on the lifetime of the pro- 
tostar, although in reality all of this mass will not be lost via the 
outflow. IHFRI looked at these limits for their outflow data and 
with our different parameters it is worth revisiting some of their 
findings. We define an outflow's mass-loss rate, M ut» as the out- 
flow mass divided by its t<j. The des truction time-scale is then, 
^des — M Qm /M 0Ut . We take M env from |Hatchell et al. |2007a| ) and 
correct for the different assumed distance. 

We find destruction time-scales in the range 0.2 to 15Myr 
for Class and 0.3 to 5 Myr for Class I protostars (see Fig. [8}. 
On average the Class time-scale is longer than the Class I: 
(f des ) = (2.5 ±0.8) Myr compared to (f des ) = (1.6±0.4) Myr. Class 
I sources have both smaller M env and M out , which counteract each 
other to yield similar £ des in both classes. There is a 67 per cent K- 
S probability that the samples from each class are drawn from the 
same population. These time- scales are smaller than HFR|s as we 



estimate larger mass-loss rates for the reasons previously discussed. 
However, they are still an order of magnitude larger than the mean 
lifetime of the embedded phase which is 1.5 to 4x 10 5 yr ([Hatchell 



|et al.|200"7a| and references therein). This leads us to HFR's con- 
clusion: outflows are not the only source of mass-loss in protostars. 
Of course accretion from the surrounding circumstellar material, 
traced by SCUBA, onto the final object may account for some ad- 
ditional mass-loss but typically the mass-loss rates via accretion 
are 10 per cent of those from outflows ( Bontemps et al. 1996). HFR 
give a thorough discussion of the potential systematic errors but two 
points are worth mentioning: (i) M env overestimates the circum- 
stellar mass remaining so we overestimate f d es- Many apparently 
single cores are in fact multiple systems. Additionally, it is well 
known that looking for clumps in two-dimensional data finds more 
massive sources than in reality exist (e.g. |Smith, C lark & Bonnell 
2008). (ii) If the destruction time-scales are truly larger than the 
core lifetimes, the problem is worse in reality. f d es was calculated 
assuming a constant mass-loss rate, whereas as we have seen it is 
likely to decrease with time, making the time- scales longer. 



5.8 Global effect of the outflows 

The outflows within every region in Perseus are energetic and ex- 
tend across significant portions of the map. Whilst it is not possible 
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Figure 7. Flux versus mass relationship corrected for luminosity with 
sources marked according to their Hatchell et al. ( 2007a) classification. 
The ordinate measures the ratio of outflow to radiative momentum flux 
(see text) and the abscissa comes from a relation in Hatchell et al. ( 2007a): 



, See also HFR 



to estimate accurately how much energy and momentum the out- 
flows actually impart to the cloud, we can estimate the maximum 
that could be shared from a simple sum of the contributions from 
each outflow. In Table [5] we list the total outflow parameters for 
the different regions. We cannot correct each individual outflow for 
inclination but a correction to the total parameters is possible, as- 
suming that the outflows are distributed uniformly, at random in- 
clinations to the line of sight. The average inclination angle is then 
(/) = 57.3 deg. For comparison, Table [SJpresents various estim ates 
of the momentum and energy from our C 18 0/ = 3— ^2 data |Pa-| 
peril. We use the average non-thermal line width of the C 18 line, 
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Figure 8. Histogram of driving source destruction time-scales. The bars 
again represent Poisson errors on the bins. 



Ont(C 18 0), as a probe of the turbulence in the bulk of the gas 
since it is substantially broadened by non-thermal motions (see |Pa-| 
perl). The turbulent energy is then E tm b ~ Mct^t and momentum 



iWh~~ Mcjnt (Wal awender, Bally & Re ipurth 2005} |Davis et al.| 
|2008| >. We estimate the mass from the C Ab O integrated intensity 
in each region, assuming LTE, T ex = 12 K and that the line is op- 
tically thin. Using this C ls O mass we calculate the gravitational 
binding energy, GM/R 2 , with the radius, R, the geometric mean of 
the major and minor axes of the cloud 'core' taken to be where the 
integrated C 18 intensity is ^1 Kkms -1 . 

Molecular outflows in NGC 1333 have sufficient momentum 
to accelerate the entire region to a few km s _1 over the cloud's life- 
time (KS00) and therefore NGC 1333 provides one of the clearest 
examples of the self-regulation of star formation. We find the total 
mass of outflowing gas is less than one per cent of the total in the 
region; the total momentum is 10 per cent of P tur b and total energy 
is around three times £ t urb but only 14 per cent of the gravitational 
binding energy. Thus there is enough energy to drive the cloud tur- 
bulence and regulate the star formation process provided the out- 
flows can be efficiently coupled to the bulk motions. Indeed there 
is a good proportion of the energy required to destroy the cloud 
in the outflows, which KS00 suggest is a likely outcome. The tur- 
bulent momentum seems too large to be provided solely by these 
outflows, although we probably underestimate the outflow momen- 
tum in CO as most of it is contained in the jet that CO traces poorly. 
Alternatively, if our momentum estimates are accurate, Davis et ah] 
( 2008) suggest that it would take a number of generations of out- 
flows to build up the turbulent momentum. 

IC348 is at the other end of the scale. Its population contains 
a larger proportion of starless cores so we can probably infer it is 
a younger region. A tiny proportion (0.03 per cent) of its mass is 
contained in outflows. Additionally, they provide less than one per 
cent of the turbulent momentum and cloud binding energy and only 
39 per cent of the turbulent energy - the lowest fractions overall. 
Only a small fraction of the turbulence in the region can therefore 
be attributed to feedback from currently active flows. 

A similar argument for LI 448 is probably misleading. Its out- 
flows clearly burst out of the natal cloud core, making it difficult to 
understand how extended high-velocity gas, outside of the cloud, 
can feedback energy and cause turbulent motions or large-scale dis- 
ruption. Nevertheless, in their entirety the outflows are massive and 
energetic, containing twice the binding energy of the cloud (also 
noted by Wolf-Chase et al. 2000), so if only a fraction of that en- 
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Table 5. Total outflow properties in each region. Estimates of the momentum, kinetic energy and momentum flux have been increased by factors of (1/ cos i) = 
2, (1/ cos 2 i) = 3 and (sin// cos 2 i) =4/71 respectively, assuming a random, isotropic distributions of outflow inclinations. 



Region 


Mass 


Momentum 


Kinetic Energy 


Momentum Flux 




(Mo) 


(M kms" 1 ) 


(J) 


(Mokms-V 1 ) 


NGC1333 


1.62 


19.4 


4.9 xlO 38 


1.8xl0" 3 


IC348 


0.05 


0.5 


l.OxlO 37 


5.5xl0" 5 


L1448 


0.75 


13.6 


5.4xl0 38 


l.lxlO" 3 


L1455 


0.10 


0.9 


1.5 xlO 37 


7.1 xlO" 5 



Table 6. Properties of the different regions from C 18 0/ = 3^2 observations. Masses are calculated assuming LTE and T ex = 12 K. The radius is the geometric 
mean of the major and minor axes of the main C ls O 'core' (see text). Turbulent momenta and energies are calculated from the average C ls O non-thermal 
linewidth, Ont(C 1s O), via P tU rb ~ Mont and £ tU rb ~ Mo^ T . 



Region 


Mass 


Radius 


Binding Energy 


ctnt(C 18 


O) 


^turb 


^turb 




(M ) 


(arcsec) 


(J) 


(kms~ 


l ) 


(Mokms" 1 ) 


(J) 


NGC1333 


439 


775 


3.5 x 10 39 


0.44 




193 


1.7 x 10 38 


IC348 


196 


440 


1.2 x 10 39 


0.26 




51 


2.6 x 10 37 


L1448 


59 


198 


2.5 x 10 38 


0.35 




21 


1.4 x 10 37 


L1455 


19 


90 


5.5 x 10 37 


0.31 




6 


3.6 x 10 36 



ergy were coupled to the cloud, the core would probably be broken 
apart. 

Finally, L1455 is similar to NGC 1333. Around 0.5 per cent of 
the total cloud mass is held in outflows, with 15 per cent of the tur- 
bulent cloud momentum and four times the turbulent energy. There 
is still not quite enough kinetic energy to overcome the binding en- 
ergy, although we have not counted the largest northwest- southeast 
flow in our calculations ( Goldsmith et al. 1984), which if included 
may well raise the total beyond the binding energy. Many of the 
outflows again stretch outside the dust and ambient emission, so 
these totals are upper limits to the energetic input into the cloud. 

Numerous previous studies have explored the total contribu- 
tion outflows might provide to the turbulent energy in other molec- 
ular clouds. In the Serpens molecular cloud the kinetic energy from 
outflowing gas is ~ 70 per cent of the cloud's turbulent energy 
and rsj 62 per cent of its gravitational binding energy ( Grave s^etaL] 
|2010| ), values perhaps most comparable to IC348 - although Ser- 
pens is closer to equipartition between its gravitational and turbu- 
lent energies. This estimate for Serpens does apply a correction for 
the outflow inclinations but not for the optical depth of the 12 CO 
/ = 3 — > 2 transition. In Section |4~T] we demonstrated this can in- 
crease the outflow mass (and kinetic energy) by a median factor 
of 3.5. Applying this correction would mean the outflows' energy 
is more than double the cloud's turbulent and gravitational bind- 
ing energy; a similar situation to that noted in LI 448. Stanke & 
Williams (2007) and Dent et al. (2009) analysed data towards two 
clustered regions forming high-mass stars: L1641-N in Orion and 
the Rosette molecular cloud respectively. In LI 641 -N a lower limit 
for the outflows' energy of ~ 30 per cent of the turbulent was found, 
again similar to our results in IC348. Once more, the outflows' ki- 
netic energy may be underestimated by over an order of magnitude 
as no correction for inclination nor optical depth was applied. Dent 
|et al.| ( [2009] ) find that their outflows (increased in mass by an aver- 
age factor to account for optical depth and inclination) have similar 
kinetic energies to their parent stellar cluster's isturb* but in total 
their energy is ~ 2 orders of magnitude smaller than E tm b for the 
entire cloud. Our results also suggest that for some clusters of star- 



forming cores outflows may be energetically significant although 
their contribution is likely to be far smaller if Perseus is considered 
in its entirety. 



6 SUMMARY 

This paper presented a survey of molecular outflows across the 
Perseus molecular cloud. We analysed large-scale 12 CO J = 3 —> 
2 datasets (over lOOOarcmin 2 ) undertaken with HARP on the 
JCMT towards four of the most active regions of star formation 
(NGC 1333, IC348, L1448 and L1455), combining previously pub- 
lished and new data. These new data across over 600 arcmin 2 of 
NGC 1333 are the largest contiguous CO J = 3—^2 maps published 
of NGC 1333 to date, extending the mapping of KS00 by a factor 
~ 10. For each detected outflow we calculate various parameters 
and use our complementary data (Paper I) to: (i) correct the outflow 
mass, momentum and energy for the optical depth of the 12 CO gas; 
and (ii) estimate the ambient cloud/driving- source velocity at each 
spatial position. Our main conclusions can be summarized: 

(i) Detection rates. Of the 65 SCUBA cores in our fields 
(Hatchell et al. 2007a), we detect outflows towards 45 (69 per cent). 
We disagree with [HFR| in outflow (non-)detection towards only one 
of our 40 common sources, suggesting that the increased sensitiv- 
ity of this survey does little better at differentiating starless and 
protostellar cores (as found by Hatchell & Dunham 2009). Indeed, 
spatial resolution is a greater limitation. 

(ii) Length and maximum velocity. We find Class sources 
drive longer and faster outflows than Class Is. (Liobe) = (140 ± 
20) arcsec and (v max ) = (17.9 ± 1.4) kms -1 compared to (88 ± 
11) arcsec and ( 12. 2 ± 1.6) kms -1 for Class and I sources respec- 
tively. This might be explained as Class I outflows have wider open- 
ing angles than Class { [Arce & Sargent 2006), so their entrained 
material is more spread out and therefore less bright. Outflows also 
become less forceful as they age (see below and|Bontemps et al.| 
1996) so their emission spans smaller velocity ranges and may be 
more easily masked by ambient emission. The ratio of the length of 



© 2010 RAS, MNRAS 000,[l]j24] 



Molecular outflows in the Perseus molecular cloud 15 



the longer to shorter outflow lobe is on average, 2.0 ±0.6. This pos- 
sibly indicates variations in the environment around the protostar. 
If, for instance, one half of the outflow breaks out of the molecu- 
lar gas closer to the protostar than the other, this outflow lobe will 
correspondingly appear shorter. 

(iii) Mass, momentum and energy evolution. Class I outflows 
have less observable mass, momentum and energy than Class 0s. 
The mass distributions for each type are not significantly different 
with means of (M out ) = (0.09 ± 0.02) M and (0.06 ± 0.03) M 
for Class and I sources respectively. The momentum and en- 
ergy distributions for the different types are significantly differ- 
ent. The Class averages are (p ut) = (0.7 ±0.2) M kms -1 and 
(£out) = (1.4±0.5) x 10 37 J compared to (0.10±0.03) M kms" 1 
and (1.0 ± 0.3) x 10 36 J for the Class Is. In the Class I popula- 
tion, SVS13 (omitted from the previous averages) is an anomalous 
source, possibly as it has a Class protobinary companion (e.g. 
|HFR] >. 

(iv) Dynamical time-scale and momentum flux evolution. We 

compute Td using the "v max " method over the whole of the out- 
flow. These are unrealistic outflow ages showing little correla- 
tion with Tboi and little difference between the two source classes: 
(r d ) = (9800 ± 1600) and (10000 ± 2000) yr for the Class and 
I average respectively. For Class I sources in particular it will 
probably underestimate the outflow age considerably, partly as we 
have trouble mapping these outflows accurately. The outflow mo- 
mentum flux computed from 1$ is nevertheless larger for Class 
sources than Class I: (Fco) = (0.8 ±0.3) x 10 -4 compared to 
(1.1 ± 0.3) x 10" 5 M kms -1 yr" 1 (again with SVS13 omitted). 
This difference is emphasized if we compute the momentum flux 
using the protostellar age derived from T^ol instead of the ob- 
served outflow dynamical time-scale: (Fq ) = (4.3 ±1.5) x 10 -4 
(Class 0s) compared to (6 ± 2) x 10" 6 M km s" 1 yr" 1 (Class Is). 
Class sources are systematically more efficient at driving out- 
flows than Class I. If F ra( j = L\, \/c is the flux expected in ra- 
diation from the central source then Fco (Class I) ~ lOO^rad an d 
Fco (Class 0) ~ 1000F ra( j. The momentum flux also correlates with 
the source mass and inversely with the source T^ \. Such a decrease 
in force between stages, with age, suggests a decline in the mass 
accretion rate (see Bontemps et al. 1996 ) or possibly a different en- 
trainment efficiency which may occur if for instance a second wind 
component has started to dominate. 

(v) Destruction time-scale. If mass is only lost from the core via 
its outflow at the current rate, then the cores will last too long. 
The average calculated destruction time-scale is a few million years 
compared to the average lifetime of the protostellar stage which is 
a few times 10 5 years. This is probably because our estimates of 
the masses of the star-forming cores are too large. 

(vi) Global effect of the outflows. We estimate the maximum 
mass, momentum and kinetic energy available to each region from 
their outflows by summing the contributions from each individual 
outflow and applying an average correction for inclination. This is 
compared to some overall cloud properties: the turbulent momen- 
tum, energy and binding energy, derived from our C 18 data. We 
find there is more energy in outflows than is observed in bulk, tur- 
bulent motions in NGC 1333, L1448 and L1455. Indeed, the out- 
flows in LI 448 are extremely powerful, containing twice the bind- 
ing energy of the region. We note the difficulty in feeding-back 
high-energy outflowing gas, often outside of the main cloud core, 
into turbulent motions. However, if it were possible with only a 
fraction of the outflow energy in LI 448, it seems likely the cloud 
would be broken apart. 
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APPENDIX A: OUTFLOW DESCRIPTION AND 
BACKGROUND 

Al NGC 1333 

A 1.1 Overview 

We present an overview of the outflows in NGC 1333 in Fig. 
|A1| Most of the high-velocity gas lies in the central regions, also 
mapped by KSOO although the full extent of these outflows is better 
traced here. The blue-shifted emission is dominated by a large-scale 
flow north- south from HRF44 and the southeastern blue- shifted 
lobe from SVS13 (HRF43). The most prominent red-shifted fea- 
ture expands into the central dust cavity, surrounded by the most 
well-know star-forming cores ( Sandell & Knee pOOT) . 
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Figure Al. High-velocity CO outflows towards NGC1333. Red- and blue-shifted CO J = 3 — >■ 2 emission has been integrated (jT^dv) from —5 to 
3kms -1 (blue) and from 12 to 18kms _1 (red) and is contoured at the levels shown (in Kkms" 1 ). The colour-scale is SCUBA 850 |n.m flux density scaled 
from to 1600 mJy beam -1 . The thick grey line marks the area surveyed for outflows with HARP. 
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A1.2 NGC1333-IRAS2 

HRF44 is better known as the infrared source NGC 1333-IRAS2, 
discovered by Jennings et al.| {1987) in IRAS data. In high- 
resolution millimetre studies, three sources are apparent: two pro- 
tostellar (2A and 2B) and one starless (2C) ( |Blake| 19961 |Looney,| 
Mundy & Welch 2000]). The protostars are also detected at cen- 
timetre wavelengths (e.g. Rodriguez, Anglada & Curiel 1999J- Two 
outflows are associated with the source, one jet- like orientated east- 
west ( |Sandell et al.|1994||Bachiller et al.|1998] > and the other shell- 
like approximately north- south ( [Liseau, Sandell & K nee 1988 ). 
KS00 point out such different structures may imply they arise 
from separate components of IRAS 2. Alternatively, a single binary 
source is possible, supported by the interferometric observations of 
Engargiola & Plambeck (1999) which pinpoint the origin of both 
flows at IRAS2A. 

A detailed picture of the outflows around IRAS2 is shown 
in Fig. |A2| The two orthogonal flows are clear. The red lobe en- 
compassing H2 knot 32, extends into the confused area around the 
central cavity and KS00 suggest it created part of this void. High- 
velocity blue-shifted gas, extends to HH13 in the south, probably 
causing a number of knots en route. The lower velocity gas is even 
more knotty and supports the continuation of the flow even fur- 
ther south ( [Davis et al.|2 008 suggest it extends 2.2 pc). However, 
HRF75 and HRF71 (Class 0) alongside HRF69 (Class I) may drive 
flows that interfere or contribute to this lobe, even though there is 
little strong outflow activity around them, perhaps just an east- west 
flow from HRF71. 



A1.3 NGC1333-IRAS4 
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IRAS4 ([Jennings et al.|1987| ) contains at least four distinct sources. 
|Sandell&K nee ( 2001) 1 resolved two components: 4A (HRF41) and 
4B (HRF42), separated by 30arcsec. Eac h of these are multiple 
objec ts unresolved with SCUBA: 4A/4A ' ( [Lay, Carlstrom & Hills' 
1995) and 4B/4B' ( |Looney et "aL]|2000| ). Another component 4C 
(HRF48), ~ 40 arcsec north-northeast of the main pair, was discov- 
ered by Rodrig uez et al.|(|1999|>. The outflow was first mapped in 
CO / = 3 -> 2 by |Blake et al.||l995) . They found a change of di- 
rection in the flow from a position angl^j(PA) of 45 deg far away 
to deg close to IRAS4A coupled with a string of symmetric fea- 
tures either side of the source, suggesting some variability in the 
ejection. In interferometric data, a small outflow from IRAS4B is 
clear (e.g. J0rgensen et al. 2007b). 

The outflows we detect are plotted in Fig. |A3| A tiny north- 
south flow from 4B is discernible along with the stronger flow from 
4A which encompasses H2 knots 32, 39 and 67. The change in its 
direction away from the driving source is clear. Beyond the high- 
velocity red lobe, Di Francesco et al. (2001) suggest the flow bends 
to the north to encompass the red-shifted gas along the eastern side 
of the entire central region (which alternatively may originate from 
a source near HH6), resulting from a variation in the jet direction 
caused by the protobinary. SiO observations by Choi ( 2005) show a 
very straight bipolar flow emanating from IRAS4A', with a weaker 
blue lobe in a north- south flow from IRAS4A. The main north- 
easterly red lobe has a dramatic change of direction to the east (just 
beyond where our CO emission ends) which Choi interprets as evi- 
dence of an obstacle in the path of the flow. Thus, it seems probable 
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Figure A2. High-velocity CO J = 3 ->• 2 outflows towards NGC 1333- 
IRAS2, overlaid on a SCUBA 850 |u.m data. Symbols mark submillimetre 
cores classified by Hatchell et al. ( 2007a) as: starless (black/white circles), 
Class (light blue diamonds) and Class I (purple triangles). Green squares 
are H2 knots ( Davis et al. 2008) and blue crosses are HH objects ( Wal awen^] 
der et al.|2005) . J T£ dv from: -10 to 1 kms" 1 (blue) and 14 to 25 km s -1 
(red). 



that our CO flow does not veer north but rather south and the ap- 
parent change of direction in the outflow near to the driving source 
is an unresolved flow from its close companion. 



AL4 The SVS13 ridge to HH12 

The ridge of dust just south of the central cavity is an exceptionally 
rich star-forming site, containing the well known SVS13 source 
( [Strom, Vrba & Strom|1976| HRF43) which is usually associated 
with NGC 1333-IRAS3 {Jennings et al.]T987] l. The HRF43 core ac- 
tually contains numerous sources, in order from north-east to south- 
west: (i) SVS13, thought to power the symmetric CO outflow (e.g. 
|Snell&Edwards|1981| > associated withHH7-ll (|Herbig|1974| >. (ii) 
SVS13B, a probable Clas s source ([Grossman et al.|1987|) and (iii ) 
H20(B), a maser source (Haschicket al. 1980| |Hirota etalT2 0Q8). 
Rodriguez et al. ( 1999} additionally found two centimetre sources, 
VLA3 and VLA20, interspersing these. 
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Right ascension (J2000) 

Figure A2 - continued Low- velocity CO gas, / T£ dv from: 1 to 5kms -1 
(blue) and 10 to 14kms _1 (red). 



Additional outflows are very confused; both KS00 and [Davis] 
et al. (2008 ) propose two further outflows from the ridge, tenta- 
tively suggested by |KS00| to be driven by SVS13B and H 2 0(B). 
The first, from SVS13B is almost north-south, its blue lobe appar- 
ent in Fig.[A3] joining the SVS13 ridge to the blue lobe of IRAS4A. 
The second from H20(B), passes close to the long blue lobe of 
IRAS2 (see Fig. \A2\ and possibly produces the chain of features 
south to HH14. There is certainly outflowing gas in the vicinity 
of HRF50 (the Class I source near the peak of the blue- shifted 
gas, which also encompasses HH7-1 1). However, HH7-1 1 are com- 
monly associated with SVS13 and not HRF50. Perhaps some of the 
nearby outflowing gas emanates from HRF50, but it is difficult and 
rather arbitrary to decompose the emission further so we do not 
calculate outflow parameters for HRF50 (i.e. in Table |Bl| ). 

KS00] propose that the horseshoe of blue-shifted gas associ- 
ated with HH12 (see Fig. |A4} is the counterpart to the SVS13B 
outflow. Its blue- shifted nature may result if the flow hits a density 
gradient and is deflected. Davis et al. (2008 ) prefer an explanation 
that infers an outflow driven by a source inside HH12 (the most 
likely candidate being HRF45, a Class I protostar). Our data do 
not contribute much further, certainly there are a lot of potential 
driving sources in the vicinity of HH12, some with outflowing gas, 




Right ascension (J2000) 




Right ascension (J2000) 



Figure A3. NGC 1333-IRAS4 as for Fig.[A2] Upper panel: High-velocity 
CO outflows, integrated from -10 to 1 kms -1 (blue) and 14 to 25kms _1 
(red). Lower panel: Low- velocity CO outflows, integrated from 1 to 
4kms _1 (blue) and 10 to 14kms _1 (red). 



although bipolar structures are not obvious. In our analysis (Section 
[5}, we attribute the horseshoe and overlapping red-shifted emission 
to nearby HRF45. 



A1.5 HH6 

The HH6 jet ( |Lightfoot & Glencross|1986[ICohen, Jones & Hereldl 
1991 ) and associated outflow ( |Liseau et al.|1988| ) to the north-east 
of SVS13 and east of the central cavity are driven by the Class 
protostar NGC 1333-IRAS7 ( [Jennings etal.|1987| HRF46). Again, 
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Figure A4. HH12 as for Fig. |A2| Upper panel: High-velocity CO outflows, 
integrated from —2 to 3kms _1 (blue) and 15 to 20kms -1 (red). Lower 
panel: Low-velocity CO outflows, integrated from 3 to 5kms _1 (blue) and 
12 to 14kms _1 (red). 
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Figure A5. HH6 as for Fig. |A2| Upper panel: High- velocity CO outflows, 
integrated from —7 to Okms -1 (blue) and 15 to 22kms _1 (red). Lower 
panel: Low-velocity CO outflows, integrated from to 5 km s -1 (blue) and 
10 to lSkms" 1 (red). 



it is not a single submillimetre core but two ([Sandell & Knee 2001), 
both protostellar (J0rgen sen et al.||2007a|). E ither one (both unre- 
solved as HRF46 in |Hatchell et al.||2007a| ) could drive the main 
flow in Fig. |A5| In addition, there is a long orthogonal flow, de- 
tected in the near IR ( |Cohen et al.|1991[|Aspin, Sandell & Russell| 
1994 ), thought to extend down the eastern side of the map from 
HH5 to the most northerly H 2 feature in Fig. |aT]( D avis et al.|2 008). 
KS00 note red-shifted CO from this flow in the south (as we do), 
but no blue- shifted to the north where they suggest there is only 
an infrared jet. We see the red-shifted CO gas continuing along the 
entire jet extent suggested b y|Davis et al.| ( [2008] ). Indeed it is this 
trail that |Di Francesco et al.| ( |2001| ) may confuse with a deflection 
north of the IRAS4A outflow. 



A 1.6 Elsewhere 

Very few of the remaining cores drive clear bipolar outflows, de- 
spite the vast majority having high- velocity gas (see Section [3}. 
One further distinct detection comes from HRF65 (see Fig. |A6| >. It 
drives a long collimated flow clear in Fig. |A1[ just to the south of 
KSOO's field. Its blue lobe crosses the southbound outflows from 
the SVS13 ridge and IRAS2 via numerous F^knots, before ter- 
minating in the H 2 knot 66 (Davis et al.l[2008] > and HH744A/B. 
The red lobe follows a more curved path, possibly extending to- 
wards HH5, although its structure in the H 2 images suggests it is a 
bow shock from the north, possibly originating from HRF46 (see 
above). Hatchell & Dunham (2009 ) also find a new outflow from 
HRF65, although their smaller map is not able to follow this highly - 
collimated outflow along such a large extent. 
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A2 IC348 

The centre of the IC348 cluster is devoid of star formation activ- 
ity with the youngest protostars concentrated in a ridge lOarcmin 
south-west (e.g. |Muench et al.|2007|>, ar ound an IR object (IC348- 
IR) first identified by | Strom et al. (|1976|). A number of dense cores 
were found in NH3 and CO by Bachiller, Guilloteau & Kahane 
( |1987| ), before [McCaughrean, Rayner & Zinnecker (1994]) discov- 
ered HH211, a molecular jet that has been the target of numerous 
interferometric studies (e.g. |GuetF & Guilloteau 1999, Chandler 



& Richer 2001). The largest existing outflow survey is by Tafalla, 
Kuma r & Bachiller (2006), in CO J = 2 1 with the IRAM 30- 
m telescope at 1 1 arcsec resolution. HFR find many more distinct 
outflows although some overlap, causing confusion. 

We present maps of the outflows in Fig. |A7| The two most 
prominent outflows are the north-south flow driven by HRF13, 
also known as IC348-MMS ( [Eisloffel et al.|2003) or IC348-SMM2 
( |Walawender eTalJ|2006l [Tafalla et al.|2006) and HH211 driven 
by HRF12. In the slow gas, the weak bipolar outflow discovered 
by |Tafalla et aT] ( [20061 from HRF15 or IC348-SMM3 can be dis- 
cerned with further one-sided flows around HRF14 and HRF101, 
which are also noted by Hatchell & Dunham (2009). 

|Eisloffel et aL| ( [2003] ) suggest the IC348-MMS outflow ex- 
tends some lOarcmin, encompassing the H2 knots in a line to its 
north along with the group of knots south, next to HH840. IC348- 
MMS drives a complex bipolar flow, in the north there is red- shifted 
gas in the same place and beyond the blue-shifted lobe. Further- 
more, slow gas in the blue lobe appears to bend around on itself. 
There are a number of different explanations for such; first the 
anomalous red emission might simply come from another outflow. 
A likely driving source is HRF17, classed as starless by |HatchelT| 
|et al.|p007a) .|HFR report a tentative outflow detection towards it, 
although they note the potential confusion with the stronger flow. 
However, the outflowing gas bends around HRF17 instead of ob- 
viously emanating from the its position. Taf alla et"aT] ([2006) sug- 
gest that this superposition of red- and blue-shifted gas could result 
if the IC348-MMS outflow lies close to the plane of the sky (e.g. 
|Cabrit, Goldsmith & Snell|1988| ) or if the outflow has broadened 
after hitting a lower density medium. The arc of blue- shifted emis- 
sion is more perplexing, given that the outflow wind must continue 
in a straight-line to encompass the northern H2 jets. Again, [Tafalla 
|et al.| ( |2006| ) explain this simply as an unrelated gas, flowing into 
a circular hole, evacuated by the cluster. In our data the blue arc 
is part of a larger structure between 4 and 6kms _1 that extends 
to the north-east, into a lower density area of the SCUBA map. 
There is some evidence of a collimated filament of molecular ma- 
terial, coincident with the southern H2 knots attributed to the same 



IC348-MMS flow by Eisloffel et al. (2003), in blue-shifted velocity 
channels between 4 and 6kms _1 and the red-shifted ones between 
11 and 13kms _1 . 

Elsewhere, there are a number of red- shifted flows with no 
blue counterpart. HRF101 has a small red- shifted flow, that lies 
in line with a number of H2 knots ( Walawender et al. 2006). The 
only outflow activity in the east resides around HRF14, whose red 
lobe encompasses HH841, an optical knot embedded in a reflection 
nebula ( Walawender et al. 2005). There is evidence of a weak blue 
lobe on the opposite side of HRF14, prominent in the 4-5 km s -1 
channel. 



A3 L1448 

L1448 harbours a rich cluster of Class protostars. IRAS found 
three sources: L1448-IRS1 to IRS3. L1448-IRS1, to the west and 
slightly north of the SCUBA cores (at 3 h 25 m 09 s , +30 °46 / 21 // ), 
has been identified as a Class I protostar (Eisloffel 2000) but con- 
tains no compact 850 \im dust emission and therefore was not ex- 
amined by |Hatchell et al.|(|2007a|).|Bally et al.| ( |2008| ) note it has a 
visual counterpart ( [Cohen & Kuhi|1979| > and optical reflection neb- 
ula association (RN013) and is thus likely to be a Class II source. 
Additionally, the c2d survey do not detect it in Spitzer's IRAcQ 
bands (J0rgensen et al. 2007a), cataloging it as 'star+dust'. How- 
ever, Davis et al. (2008) suggest this non-detection is due to satu- 
ration and associate IRS1 with a number of H2 knots and HH194. 
L1448-IRS2 is a Class protostar that |Wolf-Chase et al.|p000] > 
suggested is actually a binary (HRF30 and HRF31), driving two 
CO outflows. The third source, L1448-IRS3, consists of three or 
more Class protostars (e.g. |Q' Linger et al.|2006[ |Looney et al.| 
[2000 1 : L1448C (HRF29 also known as L1448-mm) which drives 
the best studied outflow in the region, L1448N: A,B (HRF28, a pro- 
tobinary separated by ~ 7 arcsec, Curiel et al. 1990) and L1448NW 
(HRF27). 

In Fig. [A8] we present our CO outflows. L1448C (HRF29) 
drives the highly collimated flow first discovered by Bachiller et al. 
(1990): one of the youngest and highes t velocity flows known at 



1995). A chain 



the time (out to ±701mis~^ |Bachiller et al.|l990 
of H 2 knots extends - 2arcmin south of L1448C (Eisloffel 2000). 



Towards the north, the blue-shifted flow possibly continues beyond 
the other peaks, deflected by the ammonia core containing the pro- 
tobinary L1448N:A/B ( [Curiel et al.[T999]) terminating at HH267 (to 
the north and west of our map, |Wolf-Chase et al.|2 000). The sec- 
ond outflow cuts across L1448C's and is associated with multiple 

7 Infrared Array Camera observing at 3.6, 4.5, 5.8 and 8 \xm. 
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Figure A7. CO J = 3 -> 2 outflows in IC348 as for Fig.[Al] Top: High- 
velocity CO outflows, integrated from —7.4 to 3.6 km s -1 (blue) and 13.6 
to 23.6 km s -1 (red). Middle: Low-velocity CO outflows, integrated from 
3.6 to 5.6 km s -1 (blue) and 11.6 to 13.6kms _1 (red). In the bottom panel 
squares denote H2 knots identified by Eisloffel et al. (2003) (pink) and 
IWalawender eTaLH2006) (red). 
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Figure A8. CO J = 3 ->> 2 outflows in L1448 as for Fig. |Al] Top: 
High- velocity outflows, integrated from -25 to km s -1 (blue) and 8 to 
25kms _1 (red). Middle: Low- velocity outflows, integrated from —5 to 
2kms _1 (blue) and 7 to 14kms _1 (red). 



Class objects: L1448N:A/B (HR F28) and L1448NW (HR F27). 
Previously, one (from L1 448N:B|Bachiller et al.|1990|[T^95} , two 
(from L1448N:A and B |Wolf-Chase et al.||2000| > "or three flows 
( |Eisloffel|2000l|Davis &~S mith 1996} have been suggested for this 
area, with sensitivity and angular resolution a severe limitation. Our 
maps do not add anything to this discussion; one or more of the 
sources drives an outflow with PA 140 deg. Interferometric obser- 
vations by |Kwon et al. 1(2006) find two flows: from L1448N:A at 



PA 155 deg and from L1448N:B at 1 10 deg. In channel maps there 
is a hint of two distinct flows, the red-shifted emission between 7 
and 12 km s _1 shows separate trails originating at the L1448N core. 

The IRS2 core (HRF30 and HRF31) is thought to drive 
two outflows ( Wolf-Cha se et al.||2000| ). This is implied by the 
'V -shaped blue lobe, although red counterparts are not obvious. 
HRF31 seems to drive the red lobe oriented north- south but it has 
no symmetric blue lobe. |Davis et al.| ( |2008| ) find only one flow from 
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Figure A9. CO / = 3 ->> 2 outflows in L1455 as for Fig.lAl] Upper panel: 
integration from —5 to 3 kms -1 (blue) and 9 to 16kms _1 (red). 



this core, encompassing this red lobe out to HH195. They suggest 
the doubted protostar, IRS1, also powers a roughly parallel outflow, 
reaching to HH194. Our observations indicate that IRS1 is not a 
young protostar and the outflowing gas nearby comes from flows 
powered by two Class objects in IRS 2. A plausible morphology 
is HRF31 driving an outflow at a PA of 100 deg and HRF30 one at 
130 deg. 

A4 L1455 

Outflowing CO gas was first detected in LI 45 5 by Frerking & 
Langer (1982) with the / = 1 — >> line later mapped by Gold- 
smith ett^|l984} and |Levreauitl ( |1988| ). Subsequently, |Davis et a" 
( |1997| ) mapped a small area around one of the protostars, RN015- 
FIR in CO / = 3 2 with the JCMT, while |Bally et al.| (T997] > 
examined the whole region optically. Our outflows are depicted in 
Fig. |A9| All four protostars (HRF35-39) appear to be driving out- 
flowing gas of some description. Starless HRF40 has red-shifted 
gas in its vicinity which does not have a distinct outflow shape, so 
may just be confused with other flows. 

The brightest SCUBA core is HRF35 (RN015-FIR), a Class 
I protostar identified in the IRAS point-source catalogue (IRAS 
02245+3002). It drives the clearest bipolar outflow in the map (PA 
42 deg) which Davis et al. (1997) found to have subtle variations in 



direction possibly due to the presence of a binary source. HRF36 
drives an orthogonal outflow. However, the rest of the map is com- 
plicated: there is a wide, nearly north-south outflow from HRF37, 
but most perplexing is the outflow in the north-west which seems 
to require a driving source far from the existing protostars. This 
difficulty was noted by Levreault (1988), who found no candidate 
and derived very low luminosity limits on such a source if it exists. 
Subsequent submillimetre and IR surveys have not increased the 
number of known protostars greatly; J0rgensen et al. ( 2007a| and 
|Hatchell et al.| (2007a) both found four protostars clustered around 
the three IR sources] Levreault knew. The brightest H2 emission and 
HH279 both tally with this flow and it would seem to be driven by 
a low luminosity source or is a flow from the known protostars to 
the south-east (probably from either HRF37 or HRF39). This latter 
explanation is promoted by |Davis et al.| ( [2008] ) who correlated the 
WFCAM H2 images with the outflows in the CO / = 3 — > 2 survey 
of HFR although their CO data do not extend so far north. 



APPENDIX B: SPECTRA AND PARAMETERS OF 
INDIVIDUAL OUTFLOWS 

The full version of this appendix is available as Supporting Infor- 
mation to the online version of this article. 
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Figure Bl. Spectra of outflows in NGC 1333 for sources in Hatchell et al. (2007)a. Spectra are towards the co re pe ak (black) and the peak of the integrated 
outflow emission (from vo +2 to vo + v max kms _1 ; values of vo and v max are given for each outflow in Table |B1| towards the red- shifted and blue-shifted 
outflow lobes (red and blue spectra respectively), unless distinct lobes are not apparent. Similar spectra for all the SCUBA cores examined are available as 
Supporting Information to the online version of this article. 



Table Bl. Observational parameters of the identified outflows. The full version of this table is available as Supporting Information to the online version of this 
article. 



No. 


Driving a 


Class b 




Lobes? d 


Axis PA 


Lobe Length e 


Vmax 


(kms" 




Source 




(kms" 1 ) 




E of N 
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B 


R 


B 
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HRF41 





7.11 
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7.52 
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8.24 
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HRF44 





7.62 
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111 
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131 


22 
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4B 1 


HRF44 





7.62 
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18 


120 


460 


32 
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a Driving source from |Hatchell et al.U2007a) 

Source classification (Hatchell et al. 2007ai: 0=Class and I=Class I. 
c Source C 18 0/ = 3— )>2 velocity where available or the map median where not. 
d Number of lobes: R for red- and B for blue-shifted only. 
e Maximum CO lobe length. 
f Maximum CO velocity (relative to vo). 

1 Perpendicular outflows from the IRAS2A protobinary: A approximately east-west and B north-south. 



Table B2. Derived outflow parameters, see text for details. The full version of this table is available as Supporting Information to the online version of this 
article. 



Driving M out Pout a £out b T d c F C o d L out e 

Source (10" 2 M ) (10" 1 M km s" 1 ) (10 36 J) (10 4 yr) (10" 5 M km s" 1 yr" 1 ) (10" 2 L ) 



— NGC1333 — 

HRF41 7.05 5.64 12.05 0.93 6.07 10.57 

HRF42 1.13 0.58 0.73 1.11 0.53 0.53 

HRF43 36.30 30.79 69.32 0.44 70.29 128.91 

HRF44A 9.37 5.65 9.85 0.57 9.93 14.11 

HRF44B 42.16 26.05 41.00 1.36 19.22 24.64 



a Momentum. 
b Kinetic energy. 
c Dynamical timescale. 
d Momentum flux. 

e Mechanical luminosity, i.e. £ utAd- 
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